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KINEMATICS OF THE STELLAR HALO AND THE MASS DISTRIBUTION OF THE MILKY WAY USING BHB STARS 
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ABSTRACT 

Here we present a kinematic study of the Galactic halo out to a radius of ~ 60 kpc, using 4664 blue horizontal 
branch (BHB) stars selected from the SDSS/SEGUE survey, to determine key dynamical properties. Using a 
maximum likelihood analysis, we determine the velocity dispersion profiles in spherical coordinates (ay, <jg, 
dfp) and the anisotropy profile (/?). The radial velocity dispersion profile (ay) is measured out to a galactocentric 
radius of r ~ 60 kpc, but due to the lack of proper-motion information, erg, and (3 could only be derived 
directly out to r ~ 25 kpc. From a starting value of j3 « 0.5 in the inner parts (9 < r/kpc < 12), the profile falls 
sharply in the range r « 13-18 kpc, with a minimum value of f3 = -1 .2 at r = 17 kpc, rising sharply at larger 
radius. In the outer parts, in the range 25 < r/kpc < 56, we predict the profile to be roughly constant with a 
value of f3 sa 0.5. The newly discovered kinematic anomalies are shown not to arise from halo substructures. 
We also studied the anisotropy profile of simulated stellar halos formed purely by accretion and found that they 
cannot reproduce the sharp dip seen in the data. From the Jeans equation, we compute the stellar rotation curve 
(Vcirc) of the Galaxy out to r ~ 25 kpc. The mass of the Galaxy within r < 25 kpc is determined to be 2. 1 x 10 1 1 
Mq, and with a 3-component fit to v c i rc (r), we determine the virial mass of the Milky Way dark matter halo to 
be M vil = 0.9!g;| x 10 12 M Q (R vk = 249^ kpc). 

Subject headings: galaxies: individual (Milky Way)-Galaxy: halo - stars: horizontal-branch - stars: kinematics 



1. INTRODUCTION 

Understanding the formation of stellar halos gives vi- 
tal insight i nto the formatio n history and the evolution 
of galaxies ( Maiewskil 119931: iFreeman & Bland-Hawthornl 
120081) . Under the currently favored ACDM 
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model of galaxy formation, the stellar halos are thought 
to have been built up, at least in part , by accretion 
of sa t ellite galaxies dWhite & Reesl 119781: iSearle & Zinril 
119781: iHelmi & White! 11999 1). D iscoveries of str u ctures 
like the Sagitta rius dw arf galaxy dlbata etalJ 11994 1 19951: 
iMaiewski et al.l 120031: iBelokurov et all 120061) . the Virgo 
over-densityj Juric et aLl 2008b, the Triangulum-Andromeda 
structure dRocha-Pinto et al.l 12004 IMaiewski eTal] 12004 



iMartin et al.l 2007b a nd the low-latitude Monoceros ring 



dNewberg et al.ll2002b lend further support to the idea of the 
stellar halo being formed by accretion. Other than accre- 
tion, a part of the stellar halo could also be fo rmed by in- 
situ stars. Recent hydro-dynamical simulations dAbadi et al 
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20121) of galaxy formation suggest that in the inner regions 



the stellar halo might be dominated by in-situ stars, whose 
kinematic properties are distinct from the accreted stars. 

Observational evidences of multi-component halo have 
been found in dynamical studies of SPSS calibration stars 
dCarollo et al.ll2007l 120 1 (j iBeers et al.l2012l). i n rotational be- 
havior in metallicity bins dDeason et al 11201 lab , in kinematics 
of Galactic anti-cen ter and North Galactic Po le population 
dKinman et al.l 12007). in chemical properties dde Jong et al] 
1201 Oh . and also in age d ifference between in-situ and ac- 
cr eted halo dKaliraill2012l) . There also exists a counter-claim 
by lSchonrich et al.l (1201 ll) who demonst rate that the evidence 
of retrograde signal in the outer halo in ICarollo et al.l d2010h 
sample is weak and is because of a manifestation of incor- 
rect distance estimates. Investigating, whether any signal of 
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multi-component halo could also be seen in the dispersions 
of the velocity of the halo population is important. Ulti- 
mately, studying the velocity dispersion profiles of the halos 
and comparing them with simulations might help to constrain 
the model of galaxy formation. 

Lack of proper motions and the slightly off centered po- 
sition of the Sun with respect to the galactic center poses 
a unique challenge in studying the kinematics of the stellar 
halo. At distances much larger than R Q , line of sight ve- 
locity is same as radial velocity with respect to the galactic 
center. Hence the oy profile is easy to compute at large dis- 
tances and this has been well studied; Batt aglia et al.l d20"05l) 
studied the line-of-sight velocity dispersion (<t/ os ) of a mixture 
of 240 halo objects and found that a\ os decreases monotoni- 
call y beyond r ~ 30 kpc . In the outer most halo at r ~ 100 
kpc lBattaglia et al.l d2005l) and recently. iDeason et all d2012bl) 
both see a sig nificant drop in aw valu e to ~ 50 kms" 1 . 

Conversely. iDe Propris et aild2010l) studied 666 BHB stars 
from the 2QZ Redshift Survey and found the velocity disper- 
sion profile increases at large distances. But using 910 dis- 
tant halo A-type stars, iBrown et al.l d2010l) found that there is 
a mean decline of -0.38 ± 0.1 2 kms" 1 kpc" 1 in o> over 15 
< r/kpc < 75. More recently, Kue et al] d2008l) used 2401 
BHB halo stars within 60 kpc and measured a slower decline 
in ct/oj compared to earlier studies. At small r it has been dif- 
ficult to derive the oy profile, and the o nly attempt to measure 
the o r in inner-halo was undertaken bv lSommer-Larsen et al] 
d!997l) . They find a sharp fall in oy from ~ 140 kms 1 to 
~ 100 kms" 1 at r w 12 kpc although they assume the circular 
velocity to be constant. 

In a solar neighborhood one can get useful proper motions 
and with this all th e three velo c ity disp ersions (oy, erg, o$) 
can be measured. iSmith et al] (I2009al) used the full phase 
space information of ~ 1700 halo subdwarfs from the solar 
neighborhood (< 5 kpc) and determined the velocity disper- 
sions to be (oy, erg, aj) = (143 ±2, 82 ±2, 77 ±2) kms" 1 . 
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Also, iBond et all d2010l) analyzed the proper motions of a 
large sample of main-sequence stars within the solar neigh- 
borhood (< 10 kpc) and found o> = 141 kms" 1 , ag = 75 kms" 1 , 
00 = 85 kms" 1 . A summary of the estimated values of the ve- 
locity dispersions found in the literature is given in Table Q] 

Due to the lack of reliable proper motions of the halo field 
stars, even the fundamental quantities like the tangential com- 
ponents of the velocity dispersions as well as the anisotropy 
(/3) are still badly understood beyond solar neighborhood. The 
situation is however not as hopeless as it might seem. At 
small r, using only line of sight velocity it is possible to put 
some constraints on these quantities using maximum likeli- 
hood techniques, where a model or a distribution function 
needs to be specified a priori. In an analysis of 1170 BHB 
stars ranging from 5 < d/kpc < 9 6, where d i s now th e stel- 
lar distance rather than a radius, Sir ko et aTl ( 12004b ) fit an 
ellipsoidal distributi on of velocities and find that the halo is 
isotropic. Similarly. iDeason et al.l (1201 lal hereafter Dl 1) fit a 
constant anisotropy model (power law Distribution Function) 
to 3549 BHB stars constructed from the SDSS Data Release 
7 (DR7) and find that the halo between r = 10-25 kpc is tan- 
gential whereas the distant halo within 25 < r/kpc < 50 is 
radial. Dll claims f or the tangential in ner halo are in con- 
trast with the result by lBond et alJ (120101) who found the inner 
halo to be radial in the similar region (d w 10 kpc). Dl 1 as- 
sume the distribution function (DF) to be such that the tracer 
density and the potential both are power laws. Without a pri- 
ori knowledge of the density slope their estimates will have 
some systematics. Additionally, the potential was also kept 
constant in their analysis and thus can bias the results due to 
the degeneracy between the p otential and the anisot ropy. 

In their most recent work, IDeason et al.1 (l2012al hereafter 
D12) allow both the potential and j3 to vary and thereby 
break the degeneracy, and find that the outer halo within 
16 < r/kpc < 48 is radial with (3 = 0.5+q \. Previous estimates 
of the velocity anisotropy (f3) in the solar neighborhood and 
the nearby halo are summarized in Table Q] All the above es- 
timates of an anisotropy of the distant halo is done in a large 
radial bins. Their results might be accurate for the given ra- 
dial bins and could also be the actual anisotropy of the halo 
provided the anisotropy remains nearly constant or monotonic 
through out. On the other hand, if the actual f3 of the system is 
not monotonic but has a more complex shape, then estimating 
it in the larger bin will just capture an overall property of that 
bin. 

Theor etically t here a re families of the distribution func- 
tion of iHenonl d 19731) type which result in a constant 
anisotropic system, as well with the families of models 
those have their own anisotropy profi l es, which incl ude 
Osipkov - Merri tt model JOsipkovl (119791): iMerrittl (Il985allbl)\ 
iGerhardl (fl99~fl) . iCuddefordl (1 1 99 lhJBaes & van Hesel (120071) 
and few more with the H ernquist potential-density model in 
iBaes & Deionghel ( l2002l) . The question we ask is how well do 
these anisotropy profiles match the anisotropy of the Galaxy? 
More fundamentally, how well do we know the anisotropy of 
the halo? To this end we thus compute the beta profile with 
much finer spatial resolution and without any prior assump- 
tions about density or potential. 

Another use of studying the kinematics of the stellar halo 
is to constrain the mass and the potential of the Milky Way. 
Most of the methods to estimate the mass require the knowl- 
edge of the anisotropy parameter (3. Without the unbiased 
estimate of the velocity anisotropy, constraining the mass of 



the system via the Jeans equation cou ld be uncertain by 73 % 
(for -4.5 < P < 0.44) as found by IWatkins et all 120101 in 
their estimates of mass of the Galaxy. Several other authors 
have also used this a ssump tion to estimate the mass of the 
galaxy dDehnen et alJl2006t IGnedin et alJl2010t IDeason etai] 
l2012bl) . Precaution needs to be taken while making an arbi- 
trary assumption about the anisotropy. It has been found that 
the halo stars, the satelli tes and the dark matter halo have dif - 
ferent orbital properties dAbadi et al.ll2006t ISales et alj|2007l) . 
Hence assuming a constant anisotropy f or both field stars as 
well as satellites dBattaglia et alj|2005t iDehnen et al.l 120061) 
could introduce systematic uncertainties in the mass estimate. 
Ideally, in order to break the degeneracy we must have a sep- 
arate estimation of the radial velocity dispersion, the veloc- 
ity anisotropy and underlying den sity of the population, as 
pointed out by IDehnen et al.l (120061) . 

The orbital ev o lution of the Magell anic clouds 
dLin & Lynden-Belll 119821: iBesla et all 120071) . the local 
escape speed (|Smith et all 120071) . the timing argument 
dLi & White! 120081) and ~~the study of t he kinematic of 
the tracers popula ti on (iKochanekl Il996t IXue et al.l 120081: 
IGnedin et all 120101: IWatkins et aTTl2010l) are the methods 
undertaken in order to constrain the mass of the Galaxy. 
Summarizing all these attempts to constrain the mass of the 
Galaxy, the consensus can be f ound between 0.5 - 3.5 x 10 12 
Mq. Recently using BHB stars IXue et alJ d2008l) estimate the 
mass of the Milky Way to be 0.91^ 2 g x 10 12 M Q . However, 
they make an assumption that the variation of (v c i rc /vi os ), the 
ratio of the circular to the line of sight velocity, with radius is 
same as that in simulations. In this paper, we estimate v c j rc as 
far out as possible without any assumption and then use it to 
estimate the dark matter halo mass. 

This work focuses mainly on the study of the kinematics 
of the stellar halo in order to present the unbiased estimation 
of the velocity dispersions, anisotropy parameter and circular 
velocity as a function of radius to the extent data supports. We 
use a DF which does not require any assumption to be made 
a priori about the density profile or the potential. We then 
use our measurements of velocity dispersions to estimate the 
rotation curve of the Ga laxy. The di s c and bulge mass already 
being constrained from lSofue et al.l d2009l) . we focus on con- 
straining the dark matter halo mass. Using the circular veloc- 
ity curve (v c i rc (r)) we can estimate (3{r) out to as far as a r (r) 
is available. Finally, we compare our results with simulations 
in which the halo is formed purely by accretion. 

This paper is organized as follows: in Section [2] we dis- 
cuss the theoretical aspect of our analysis, the methodology 
adopted, and the details about the sample; in Section [3] we 
present our result for the velocity dispersions, confer the re- 
sults between the alternative models, and investigate the con- 
tribution of the halo substructures ; in Section [5] we present 
our estimation of the mass of the Galaxy. Results are then 
compared with the simulations in Section [6] In Section [7] we 
present our conclusion and discuss our result. 

2. THEORY AND METHOD 

We are interested in calculating the velocity dispersions 
(<T r ,CTg,(T0) for the stellar halo. However, the data we have 
is line of sight velocities. To proceed we need to make some 
assumptions about the position and the velocity of the Sun 
with respect to the galactic center. We assume R Q = -8.5 kpc, 
the velocity of the local standard of rest (LSR), Vlsr, is taken 
to be IAU adopted value = 220 kms" 1 , and the solar motion 
with respect to LSR (U, V, W) Q = +11.1, +12.24, +7.25 in 
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Table 1 

Velocity Dispersions and Anisotropy Reference Table 



Sample (number) 


distance (kpc) 


o>,o- e ,o> (kms" 1 ) 


anisotropy (/3) 


Reference 


BHB (~ 700) 


(r < 20, r > 45) 


o> = (140,90- 110) 


(radial(0.5),tangential(- 1 .3)) 


Sommer-Larsen et al. 1997 


BHB (1170) 


5 < d < 96 


101.4±2.8, 97.7±16.4, 107.4±16.6 


nearly isotropic 


Sirko et al. 2004b 


BHB (1933) 


16<r<48 




radial (0.5^ ? 8 ) 


Deason et al. 2012a 


Subdwarfs (~ 1700) 


d < 5 


143 ±2, 82 ±2, 77 ±2 


radial (~ 0.69) 


Smith et al. 2009b 


MS (10 s ) 


d~ 10 


141,75,85 


radial (~ 0.67) 


Bondet al. 2010 


BHB (3549) 


(10<r< 25,25 <r<50) 




(tangential (-0 . 6) ,radial( 0.5)) 


Deason et al. 201 la 



kms" 1 ( ISchonrich et a l. 2010). Spherical and heliocentric co- 
ordinate system are expressed in terms of (r, 6, (j>), and (d, I, 
b) respectively. 

2.1. Distribution Function 

The distribution function (DF), f, is defined such that f(x,v) 
d 3 x d 3 v is the probability of finding a randomly picked star 
in a phase-space volume d 3 x d 3 v. In general, we consider the 
stellar halo as an anisotropic spherical system. The anisotropy 
is defined as 



/? = !■ 



2a}. 



(1) 



ov,<7e and (70 being the velocity dispersions in spherical co- 
ordinates, and it describes the orbital structure of the system. 
The values of this parameter range from -oo for purely circu- 
lar trajectories to 1 for purely radial orbits. 

Families of DFs that generate the collisionless anisotropic 
spherical systems with con stant or varying veloc i ty anis otropy 
can be found in detail in iBinnev & Tremain e (2008). One 
such distribution function with constant anisotropy is 

f(E,L)=f(E)L- 2 ?. (2) 

Here E = $(r)-(v 2 /2) is the relative energy per unit mass and 
L is the modulus of the angular momentum vector per unit 
mass. Re cently, a DF given b y Equation (|2]i with the energy 
term from IE vans et al.l(l 19971) : 

f(E,L) OC £(«7-2)/7)+W7)-3/2 L -2,3 ; (3) 

was used bv lDeason et al.l d201 lab and [Deason et aTj d2012al) 
to study the rotation, anisotropy and mass of the Galactic halo. 
The parameters a and 7 are the logarithmic slopes of the den- 
sity (p oc r~ a ) and potential ($ oc r" 7 ) respectively. Hereafter, 
we refer to this function as Dl 1 DF. 

If one is interested in deriving the dispersion profiles, a 
simple distribution function that one can use is the Gaus- 
sian velocity ellipsoidal distribution function (GVE DF). The 
GVE DF has b e en use d in the context of the stellar halo by 
iFrenk & White] (fl98(ih usin g globular clusters as tracers and 
by ISirko et alT(l2004bl) and I Smith etail (l2009bl) using halo 
stars. A GVE DF with rotation about z-axis is given by 



/(r,v): 



P(r) 



exp 




11) 



(4) 

The DF as given by Equation (|4]i assumes that the veloc- 
ity ellipsoid is perfectly aligned with the spherical coordi- 
nates, but in general the velocity ellipsoid ca n hav e a tilt. Us- 
ing ha lo subdwarf stars ISmith et alT (12009 al) and Bond et al. 
(2010) have found that the tilt is small and consistent with 
zero. Hence, it is safe to ignore the tilt while computing ve- 
locity dispersions. 



2.2. Parameter estimation 

The proper motion information of the stars in the stellar 
halo beyond solar neighborhood (r > 10 kpc) is not accurate 
enough to properly constrain the tangential motions. Never- 
theless, our position in the Galaxy still makes it possible to 
constrain these quantities by utilizing the tangential informa- 
tion carried by the line-of-sight velocities of the stars. How- 
ever, for that we need to marginalize the distribution function 
over the unknown quantities, which in this case are the tan- 
gential components (v/, v;,). The expression for the marginal- 
ized DF is given by, 



F(l,b,d,vi os \a n ae, 



t)- 



fdvidv b . 



(5) 



We use maximum likelihood method to estimate the model 
parameters. The log-likelihood function which we maximize 
is given by 

n 

£(l,b,d, vi os 1 0> , a g , 0-0 , v rot ) = logf\/ ; , , , vi OSi ) , (6) 

where n is the number of stars in the system under study. We 
use Markov Chain Monte-Carlo (MCMC) with the Metropo- 
lis Hasting algorithm (MHA) to obtain the posterior distribu- 
tion. We quote the central values of the velocity dispersions 
(ov,<7e,(T0) as our initial estimates and 16 and 84 percentiles 
as the error associated. Note, for the GVE DF the density term 
p(r) in Equation (@]i is not a function of model parameters and 
hence it does not have any effect on the likelihood analysis. 

Once the radial velocity dispersion o r a nd the anisot ropy 
parameter /? are evaluated, Jeans equation jJeansl[T9T5b can 
be used to estimate the circular velocity v c i,- c of the spherical 
system in equilibrium using the relation 



dlnp diner, 2 ^ 
dlnr dlnr 



(7) 



where p oc r~ a is the density of the tracer population, which 
implies d In p/d In r = -a. Through out the analysis we assume 
the density to be double power law with a = 2.4 (r ^ 27 kpc) 
an d a = 4.5 (r > 27 kpc ) in agreement with the re cent works 
by lDeason etail (1201 lbl) and lWatkins eTail J2009I) . 

For systems with constant anisotropy and a given v c ; rc , the 
solution to the differential equation (0 subject to the bound- 
ary condition lim^oo aj = reads 



1 



dr'r' 2l3 p(r')(d<f>/dr') (8) 



Assuming density and anisotropy are known we can use this 
solution to estimate a r as a function of r. 

2.3. DATA: BHB stars 
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Figure 1. Radial distribution of BHB stars in galactocentric coordinates. The 
distribution has a peak at around 16 kpc. Most of the stars are found to lie in 
the range 10 < r/kpc < 25. 

Being luminous and having nearly constant magnitude 
BHB stars are ideal for studying the stellar halo, and this is 
what we use in our study. We use the BHB catalog published 
by XI 1 for our analysis. The catalog comprises of 4985 BHB 
stars obtai ned from Sloan D igital Sky Survey (SDSS) Data 
Release 8 (lAihara et al.112.01 ll) . The stars were selected by im- 
posing limits on color and Balmer line profile measurements. 
Imposing limits on Balmer line profile measurements allow 
one to remove the main sequence stars and Blue Stragglers. 
Further details o n BHB candidate selection can be found in 
IXue et all d2008l) and references therein. To avoid contamina- 
tion from the disk stars, we restrict our analysis to stars having 
a distance |z| > 4 kpc from the galactic mid-plane. As men- 
tioned earlier, the way our likelihood function (Equation [5]l 
is laid out, this cut in distance above the plane will not intro- 
duce any bias. No velocity limits have been imposed to obtain 
the sample and thus for the purpose of kinematic studies the 
population of BHB stars we select can be considered to be 
kinematically unbiased. 

For the stars that we study the angular position is known 
very accurately, but the distance and radial velocity have some 
uncertainty associated with them. To get more accurate dis- 
tances we recalibrate XI 1 distances using a col or-magnitude 
relation derived for the same population from iDeason et alj 
(1201 lbl) . The estimated dispersion in g-band magnitudes is 
0.13, eq uivalent to a distan ce uncertainty of 6%. For the 
SEGUE (lYanny et alJ 120091) radial velocity measurements, 
94% of our sample have an uncertainty of less than 8 kms . 

The galactocentric radial distribution of the final BHB sam- 
ples is shown in Figure Q] It can be seen that the distribution 
peaks at around 16 kpc. Most of the stars are found to lie in 
the range 10 < r/kpc < 25. 

3. VELOCITY DISPERSION PROFILE OF THE HALO 

We study the kinematics of the halo in radial bins to obtain 
radial profile of the model parameters (a ri ue,a^) and also 
(3. Using only line-of-sight velocity information the tangen- 
tial components, <jq and a^, are difficult to constrain except 
in the very inner regions of the halo. However, a r can be 
well constrained both in the inner and the outer halo. This 
means that relatively larger number of stars (>1000) per bin 
are required to estimate <jg and o$ as compared to a r . Given 
that we only have about 4000 stars this means that we cannot 



Figure 2. Radial ve locity dispersion in radial bins. The black dashed line is 
the cti os profile from Xue et al. 12008). Black dashed-dotted and red dashed 
lines are the Sommer- Larsen profiles given by Equ ation (5) for the fitting 
parameters taken from Sommer-Larsen et al. 1 1997) and from the fit to our 
estimated values of o> respectively. 

measure the ae and profiles with sufficient spatial reso- 
lution. Hence, we employ two different binning schemes or 
estimators, one for radial velocity and the other for tangen- 
tial velocity. The estimators are; the equi-populated estima- 
tor (hereafter EPE) and central moving estimator (CME). In 
EPE the data is binned radially with each bin containing equal 
number of particles and this is used for computing o>(r). In 
CME a set of equi-spaced positions in r are chosen and then 
at each position an equal number of points either side of the 
chosen central value are used to estimate the desired quantity. 
We use the CME for computing ae(r) and a<p{r). The crucial 
difference between the two schemes is that while the bins are 
non-overlapping in the former, in the latter they can be over- 
lapping. In EPE the spacing between the bins is directly pro- 
portional to the number of particles in each bin «bin- Hence, if 
the desired quantity can be estimated with sufficient accuracy 
employing small «bin, then EPE is the desired method. How- 
ever, if this is not the case then it is better to use the CME 
method as the spatial resolution is not directly dependent on 

"bin El 

Finally, for our data the number density of points in r is 
highly non-uniform, and hence it is not accurate to assume 
that the desired quantity has been estimated at the center of the 
bin. To alleviate this number density bias, for both schemes 
we compute the final position of the bin as the mean r of the 
points in the bin. 

3.1. Radial velocity dispersion profile ( a r { r) ) 

Here we focus on the nature of the a, (r) profile of the halo 
for which we adopt the EPE method with «bin = 400. As ex- 
plained previously, o> can be measured out to the extent of 
the data (r ~ 60 kpc). The values of o> obtained from the 
likelihood analysis are given in Figure [2] and the error bars 
represent the 1 a confidence interval and are determined from 
the likelihood fitting. 

We find that the radial velocity dispersion, o> at the Sun's 
position (R0 = 8.5) is 145.6 kms" 1 . However, beyond the 
solar neighborhood a r sharply decreases until r ~ 15 kpc, af- 
ter which it decreases much slowly and approaches a value 

1 For the effect of the bin size, see appendix lAl 
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of around ~400 kms" 1 at 56 kpc. The error bar in a r for 
the inner halo population is large and is mainly because the 
vios contains less radial velocity information as compared to 
outer parts. The overlaid black dashed line is the linear ap- 
proxim ation (w 111 — 0.31r) for <r los profile from Kue et al.l 
(120081) . Additionally, other previous at tempts to fit the pro- 
file for ay in t he ou t er parts of the halo (Battaglia et al 1 120051: 
iGnedin et al.l 120101: iBrown et al.l 120101) also found profiles 
similar to IXue et aL | <|2008j), with slightly varying slope and 
normalization. All these profiles are reasonable estimates of 
a r for the outer halo (d ^Rq) where the assumption o> w <7] os 
holds. In the inner halo (r < 15 kpc) however the approxima- 
tion breaks down and ay strongly deviates from <7i os . It can 
be seen that the deviation of oy from o\ os increases as one ap- 
proaches the center and at R© the deviation is as high as ~ 40 

k ms" 1 . 

iSommer-Larsen et al.l 0997) provide a functional form for 
fitting the ay profiles which is given by 



al- 



7T 



_7r -1 


(?)] 


tan 1 


2 





(9) 



This has a shape which is similar to our oy(r) profile and we 
fit it to find (T = 94.5 kms" 1 , a+ = 122.3 kms" 1 , r = 13.2 kpc, 
and I = 2.6 kpc. In this function the fit parameter <to gives the 
asymptotic value that ay achieves in the outer halo, whereas 
(i7o + c+) gives the approximate value for oy in the inner 
halo. The fit parameters ro and I determine the turn-off point 
and the steepness of the transition of the profile respectively. 
A lower value of I gives a steeper transition, this can be seen 
from the comparison between the Sommer-Larsen fit (I = 7.5 
kpc) and the red line in Figure|2]which is our fit having smaller 

3.2. Tangential velocity dispersion and anisotropy profiles 
(ae(r), a^(r) and fi(r)) 

There have been few attempts to constrain the tangential ve- 
locity dispersions and most of the studies are either restricted 
to solar neighborhood (ISmith et al.ll2009bt iBond et al.ll2010l) 
or for the overall system (Si rko et al.ll2004al) . The tangential 
velocity dispersions not only provide the information of the 
anisotropy of the stellar population (through Equation [U but 
together with ay also helps to measure the mass distribution 
of the Galaxy. 

We estimate the tangential velocity dispersions (ag and a^) 
using CME with «bin = 1200 stars. Given the quality of the 
data, we are only able to measure ag and er^ out to ~ 25 kpc. 
Our estimates of ag and a^ are shown in Figure |3]5,c by the 
black dots with error bars. For uniformity, we also estimate 
the a r with this binning scheme and this is shown (only out to 
~ 25 kpc) in Figure^ by the black dots with error bars. 

In general the tangential components ag and a^ near the 
solar neighborhood are found to be comparatively lower {ag = 
85^ kms" 1 , = 95^g kms" 1 ) than the radial dispersion ay. It 
can be seen in the Figure [3j? and c that there is a sharp rise 
in the values of ag and ay/, at r = 17 kpc. Beyond this 
falls whereas og rises, given the large uncertainties and low 
number of independent bins it is unclear if this is a real or a 
spurious trend. 

By substituting the estimates of the tangential and the ra- 
dial velocity dispersions obtained using CME («bin = 1200) 
from the above analysis into Equation <£TJ we compute the 
corresponding values of the anisotropy constant in the re- 
spective bins. As shown in Figure [3}l, the halo within 12 
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Figure 3. Velocity Dispersions and Anisotropy profiles: From top to bot- 
tom is oy,<79,<70, and profiles of the stellar halo estimated in the radial 
bins. The diamond and the round markers are the results for the two binning 
schemes , namely the EPE and the CME respectively. Note that the last ra- 
dial bin marked with the open diamond contains the remaining stars. The 
diamond markers in plot (a) is just shown for the ease of comparison. 
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kpc has ~ 0.5 whereas the outer halo beyond the turn-off 
point is nearly isotropic. We discover a significant drop in 
the profile at r = 17 kpc. Here the halo is strongly tangen- 
tial with = —1.2. We later confirm that the trend observed 
in anisotropy is neither due to the manifest ation of the sys- 
tematic s introduced by the chosen m odel ( 33.31 ) nor due to 
presence of the halo substructures ( §I3.41 >. It is also shown 
in the appendix [B] that assumed Vlsr and R have negligible 
effects upon these estimates. The probable reasons for this 
sudden turn-over in the properties of the halo are discussed in 
the conclusion. 

We know that the consecutive CME bins overlap in radius 
and thus the dispersion profiles demonstrated in Figure [3] is 
a smoothed version of the actual dispersion profiles of the 
halo. However, to check for any systematic associated with 
the choice of the binning scheme, we also estimate a r ,ag, cr^, 
and in traditional equi-populated (EPE) bins («bi n = 700). 
The measured values in these bins are shown by the diamond 
points in Figure [3] If the number of stars per bin is less than 
700 it is difficult to constrain ag and o$. Even with «bi n = 700 
we were able to constrain the tangential motion only till 16 
kpc (first three diamond points). Hence, we construct the 
last bin by grouping all the stars beyond 16 kpc into one bin 
(rightmost diamond point). More importantly, with this bin- 
ning scheme we are only interested to see whether we ob- 
tain the corresponding dip or rise (depending on the param- 
eter of interest) seen in Figure [3] or not. We find that except 
for the right-most diamond points all of other diamond points 
in the figure are in agreement with our previous estimates of 
the dispersions (given by black dots). The rightmost diamond 
points are calculated in a huge bin with more than 50% of 
the total sample. Particularly for a$, a^ and given the non- 
monotonic trend they have, hence we do not claim the last 
diamond point is the correct estimate of anisotropy at r ~ 35 
kpc. 

None of the uncertainties quoted in the above estimates of 
a r , ag and a^ include the uncertainties in distances and ra- 
dial velocities. As mentioned earlier, the errors in distance 
and radial velocity are quite small, and convolving the model 
(Equation 0|i with the error functions should not change the 
results. 

3.3. 0(r) from fitting Distribution Function (f(E,L)) 

It would be interesting to see whether the profile pre- 
sented above, in particular the dip seen at r = 17 kpc, is an 
artifact of our chosen model (GVE) or a real inherent feature 
of the Galactic halo. In order to pursue it, here we explore 
the effect of the chosen model on the determination of the 
velocity anisotropy (f3). For the comparative study, the al- 
ternative model we choose is the Dll DF (Equation [3}. We 
consider anisotropy (0) and potential (= ^o?" 7 ) as free model 
parameters and constrain them using the maximum likelihood 
method. 

First, we compare the theoretical properties of the DFs at 
our disposal namely, GVE and Dll DF. In Figure |4] we plot 
the theoretical LOSVDs of models along two separate line- 
of-sights. We also show the LOSVDs in Figure|4]for two dif- 
ferent distances representing the inner halo (d = 15 kpc, given 
by red lines) and the outer halo (d = 50 kpc, given by blue 
lines). For the inner-halo we assume a = 2.4 and assign radi- 
ally biased anisotropy (0 = 0.4) whereas for the outer-halo we 
assume a = 4.5 and as sign tangential l y bias ed a nisotropy (0 = 
-2.5) in accordance to lWatkins et al.l (120091) and lDeason et al.l 
d201 lbl) estimates for a. Note, the density normalization at 



the break radius (r = 27 kpc) is assumed to be equal. For the 
assumed constant potential the solid lines in Figure [4] (both 
left, right panels) are the LOSVDs obtained by adopting the 
Dl 1 model and dashed lines are the LOSVDs of our GVE 
model. Recall that our model does not take directly but de- 
mands the information of the velocity dispersion components 
(ay, ag , a^) individually. To make the LOSVDs obtained from 
both the models comparable we estimate a r ,ag,a,p from the 
set of values of 0,a,& o and 7 chosen to obtain LOSVDs of 
Dll DF. For an assumed potential power law we use Equation 
(0 to first calculate ay, to put in our model. For an assumed 
0, substituting this a r in Equation (Q~|) gives the corresponding 

value for a t (= Jo^ + a^). It can be seen in the figure that 

for all the four cases LOSVDs obtained from both the mod- 
els match well. Naively, from these perfect matches of the 
LOSVDs at different line-of-sights and distances it can be an- 
ticipated that the estimation of 0(f) with both models should 
also match. 

Now we estimate the 0(f) using ex actly the same sample of 
BHB stars in same radial bins as in 33.2K CME. «bin = 1200) 
but with Dll DF. Blue points in Figure [5] demonstrate the 
profile obtained by fitting Dll DF. Here we use brute-force 
grid based analysis to constrain the model parameters 0, $ 
and 7. To facilitate the comparison our estimates from Figure 
[3]l are over plotted in Figure [5] and are shown by the black 
points. From Figure[5]it can be seen that within the range of 
uncertainties the measured values of with both models (Dll 
DF and GVE model) agree. However, a slight bias can be 
seen in the sense that obtained from GVE DF is in general 
higher. The reason for this discrepancy lies in the fact that 
obtained from the Dll DF has a dependence on a. Hence, 
unless the underlying a value of the sample is exactly known, 
a mismatch is expected. The estimated value of increases 
with the adopted value of a (see Fig-3 Dll). This suggests 
that the actual value of a is even higher than the one that is 
adopted here (2.4). 

In order to give an estimate of the quality of constraints ob- 
tained from the likelihood analysis we display the likelihood 
contours in the parameter space in the top left and bottom two 
plots in Figure [6] for a bin centered at r = 16.93 kpc where 
the maximum dip in the profile was seen in Figure[5] Addi- 
tionally, in the corresponding bin, the top right plot in figure 
demonstrates the posterior distributions of the model param- 
eters ay, erg and a^ of the GVE model obtained from 5 x 10 5 
MCMC random walks. It can be seen that even at a distance 
of just twice of Rq, ag and a^ distributions are quite broad 
as compared to ay distribution; this is the reason for the large 
uncertainty in the value of as we move outwards in r. 

3.4. Effect of the halo substructures 

There is now enough observational evidence to suggest that 
the stellar halo is highly structured partic ularly as one moves 
outwards into the halo dBell et al] 120081) . Usi ng clustering 
algorith ms on simulated N-body stellar halos, ISharma et alj 
d201 lbl) find that the fraction of material in substructures in- 
creases monotonically as a function of distance from the cen- 
ter and at around 65 kpc can be as high 50%. Hence, while 
studying the kinematic properties of the halo should we in- 
clude the substructures or exclude them? If the kinematic 
properties of a sample are dominated by a few massive ac- 
cretion events then one should exclude them. However, in- 
spite of being highly structured if the sample is a superposi- 
tion of large number of events with none of them being in- 
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Figure 4. Comparison of the marginalized DF given by Equation (3} obtained for Dl 1 models and GVE distribution along two different line-of-sights. Solid 
line is the LOSVD for the D 1 1 DF whereas dashed line is the LOSVD for the GVE model. Red and blue lines represent LOS VD for two different distances (d) 
15 kpc and 50 kpc respectively. The potential is assumed to be a power law given by <J>or~ 7 , where <I>o (in knrs ) is the potential normalization. Potential slope 
(7) is assumed constant and equal to 0.35 for all the cases. Left: LOSVD along (l,b) = (50°, 50°) Right: LOSVD along (l,b) = (230°, 50°). 




Figure 5. Anisotropy (ft) estimates in the CME using D 1 1 model (blue 
points) and GVE model (black points). Each bin consists of 1200 stars. 
Anisotropy estimates with GVE distribution is done with MCMC technique 
whereas estimates with D 1 1 model is done with the brute-force grid based 
analysis. Assigned asymmetric uncertainties are la confidence intervals ob- 
tained from likelihood fitting. 

dividual ly too dominant then it is best to include them. Re- 
sults of iSharma et al.l d201 lbl) on simulated halos show that 
for the range of radii that spans our BHB stars (r < 40 kpc) the 
amount of material in substructure should be less than 20%. 
So we expect the substructures to have a minor effect on the 
kinematic properties that we have derived. However, it is still 
important to check if this is true. 

To study the effect of substructures on our estimation of the 
dispersions we mask two prominent features of the halo that 
contaminate our sample, namely, the Sagittarius stellar stream 
and the Virgo over-density. Cuts we impose include the Lam- 
bert Equal-Area projection cut as given in iBell et all (120081) 
and an additio nal cut in equa t orial co ordinates (RA and DEC) 
suggested by iDeason et alJ (1201 lbl) . We mask the region 
within < X (abscissa of the equal-area projection)< 30, 
where X is given by 63.63961 [2(l-sinfo)]V2; and 0° < RA < 
50°, and -30° < DEC < 0°, which is purely a geometric cut. 
These stringent cuts reduce our final sample to 2975 stars. A 
proper phase-space masking of these structures will be revised 
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Figure 6. The posterior distributions of the parameters for bin centered at 
r = 16.93 kpc. Upper Right: The posterior distributions of the GVE model 
parameters, o>, ag and cr^, obtained with 5 X 10 5 MCMC random walks. In 
the inset is the derived distribution of the /3 parameter. Upper Left and Lower 
panels: The joint likelihood contours of the Dll model parameters /3, poten- 
tial normalization ($o) and potential slope (7) are obtained with brute force 
analysis. The outer contour displays 1 a region whereas the inner contour 
demonstrates a region of 50% confidence interval. Cross hair corresponds to 
the point where the likelihood is maximum. 

in the future work. 

In Figure|7]we present our result obtained after masking the 
substructures. As masking reduces the sample size almost by 
half, we employ CME with «bin = 500 stars only, instead of 
«bin = 1200 as was done earlier with unmasked data, to avoid 
excessive smoothing of the estimated profiles. Figure|7]shows 
that the masking of substructures have little effect on the esti- 
mation of velocity dispersion profiles, the f3 profile is almost 
unchanged. This alleviates the concern that perhaps the turn- 
over points in the velocity dispersion profiles and the dip in f3 
profile, seen in Figure [3] in the region r = 13-18 kpc could 
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Figure 7. Effect of the halo substructures, namely, the Sagittarius stellar 
stream and the Virgo over-density. Each CME bins contain 500 stars. Black 
and red points are our result for the masked and the unmasked halo respec- 
tively. Error bars quoted are lcr credibility interval. 

be due to the dominance of halo substructures. 

4. COMPARISON OF ANISOTROPY ESTIMATES IN Dl 1 AND D12 
RADIAL BINS 

In their recent work Dl 1 and D 12 fit a distribution function 
of the form given by Equation (0, to the BHB samples ob- 
tained from SEGUE survey in order to estimate the model pa- 
rameters. The models adopted are constant anisotropy models 
given by Equation In Dl 1, the potential is assumed to be a 
power law (oc r~ 7 ) and with a constant index (7 = 0.35). Later 
in D12, they break the degeneracy present in their model and 
consider the potential normalization ($0), potential slope (7) 
and anisotropy ((3) as free parameters. Note, in Dl 1 there is an 
additional parameter specifying rotation (odd part of the DF) 
but this was dropped in D12, as they were not focusing on 
rotation. The methodology applied to measure the model pa- 
rameters is similar to ours which involves marginalizing the 
DF over the tangential velocities to derive the line-of-sight 
velocity distribution (LOSVD); fitting the LOSVD to the data 
using the maximum likelihood method and in return obtaining 
the best estimates of the model parameters. 

4.1. Anisotropics of inner and outer halo by Dll 

Dl 1 study the rotation and the anis otropy of the BHB sam - 
ples taken from SDSS Data Release 7 (lAbazaiian et al.l2009t) . 
First, in order to construct the sample used by them we query 
SDSS DR7 database to select the candidate BHB stars us- 
ing the color and the stellar parameters ranges given in D 1 1 . 
Like them we also mask the Sagittarius dwarf galaxy which 
reduces the original sample size by 40% to ~ 3500. 

Dll measure the anisotropy of the halo in radial and metal - 
licity bins. In the inner halo (10 < r/kpc < 25) both metal- 
rich ([Fe/H]> -2) and metal poor ([Fe/H]< -2) stars are 
found to be tangential with (3 ~ — 1.2 and ~ -0.2 respec- 
tively. In the corresponding metallicity bins, the outer halo 
(25 < r/kpc < 50) is found to be radial with /3 of ~ 0.4 and 
~ 0.5 respectively. Since they do not give the estimates of (3 
in combined metallicity bins for the inner and outer halo, we 
estimate them here using the same methodology as adopted 
by them. For the inner halo we find (3 = -0.62 (tangential) 
and for the outer halo we find (3 = 0.41 (radial). These esti- 
mates are consistent with Dl 1 results, if we combine their low 
and high metallicity (3 values in each radial bin. 
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Figure 8. Posterior distribution of velocity dispersions for D12 data set 
within 16 < r/kpc < 48 using a GVE model. The value of v^sr here is 
taken to be 240 kms~' to keep it same as in D12. 

The inner halo within the solar vici nity has been found t o be 
radial in studies of halo subdwarfs bv lSmifh et al.l ( 2009b) and 
in studies of 10 5 main sequence stars by iBond et al.l (120101) 
(see also Figure [3}I). Hence, a tangential inner-halo as pre- 
dicted by Dl 1 is surprising. The first thing to check is if the 
Dll result is due to some of the assumptions made by them. 
For example, in Dll the logarithmic density slope was as- 
sumed t o be constant and equal to -3.5. Later on lDeason et alj 
d201 lbl) conducted a detailed analysis of the BHB stars to es- 
timate their density profile and found that the profile is of the 
form of a broken power law, the inner-halo (< 27 kpc) hav- 
ing a profile index of -2.3 and the outer halo having a profile 
index of -4.6. This i s in good agreement with the findings of 
IWatkins etldl (120091) that the halo within 25 kpc has the pro- 
file index -2.4 and beyond which it is steeper with an index 
-4.5. 

If we update the density profile index in the Dll case with 
the above values then we expect the inner halo which is al- 
ready tangential to become even more tangential and the outer 
halo which is already radial to become even more radial. This 
is because j3 has a dependence on the adopted value of a as 
shown by Dll (Figure 5). In general f3 increases with an 
increase in the adopted value of a. Another effect that can 
potentially bias the results is the fact that the potential param- 
eters ($0,7) have been kept fixed in the Dll analysis, i.e., 
the degeneracy between (3 and potential has not been broken. 
After analyzing data in finer bins and breaking the degeneracy 
among the model param eters we do see that in Figure[5]the in- 
ner halo is radial (from ^3.3b as was also found with the GVE 
model. Hence the assumption of a fixed potential can bias the 
estimation of (3. However, a much more apparent reason for 
the discrepancy is as follows. Clearly, from Figure [3] the ra- 
dial bin from 10-25 kpc will encompass the stars within 13-17 
kpc which are predominantly tangential. Since the probabil- 
ity density of stars in radius also peaks at around 16 kpc (see 
Figure[TJ we anticipate the overall (3 to be tangential. To con- 
clude, the tangential behavior of the inner-halo seen by Dll 
is most likely due to the large radial bin size adopted by them. 

4.2. Anisotropy at 16 < r/kpc < 48 seen by D12 

D12 re-calibrate the distances of BHB stars using the 
color-magnitude relation given by iDeason et al.l (1201 lbl) and 
then selec t stars within 16 < r/kpc < 48 ( 1933 stars) from 
IXue et all (120111) BHB samples. They fit the Dll model to 
study the nature of (3 in the outer halo. Unlike Dll, as men- 
tioned previously here they break the degeneracy present in 
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Figure 9. Maximum likelihood analysis of the anisotropy parameter, j3, for 
stars in the radial bin 35 < r/kpc < 84. Solid black lines are the 1<t confi- 
dence region. Density and potential slopes are taken to be 4.5 and 0.35. 



their model and consider the (3, <I> and 7 as free parameters. 
Thus while they fit the model they simultaneously estimates 
all three parameters. They find j3 = 0.4^2 for a = 4.6. Using a 
model allowing for oblateness (q = 0.59) they find j3 = 0.5+q \. 
If we apply the GVE distribution function to stars in the range 
16 < r/kpc < 48 we find /3 = -0.14^;^ (o> ■■ 



: 97.3^; kms" 



<rg = 122.7*330 kms -1 and = 78.5^07 kms -1 ). This more 
or less looks like the mean value of (3 in this range (see Fig- 
ure |5J provided we take into account the fact that the num- 
ber density of stars peak at around r = 16 kpc. Although the 
D12 value is still within our 1 a region, our predicted value is 
lower than D12. It can be seen from Figure[5]that j3 is not con- 
stant in the range 16 < r/kpc < 23. It increases from being 
tangenti al to isotropic. Using the Dll DF also gives similar 
result ( 93.31 1, Beyond this range we cannot directly measure 
(3, but by deriving a best fit circular velocity profile and mak- 
ing use of o> profile which is available till r = 56 kpc we can 
predict f3, and this is returned to in section [5] However, an 
assumption about a also has to be made. Beyond, r > 27 
kpc the density slope has been shown to change from -2.4 to 
-4.5. Adopting a steeper density slope increases the value of 
(3. For a = 4.5 we find (3 ~ 0.5 for r > 27 kpc; this is more in 
agreement with D12. To conclude, the D12 value of j3 = 0.4 
for 16 < r/kpc < 48, although derived for a sample which is 
dominated by stars within r < 27 kpc, is not appropriate for 
the range 18 < r/kpc < 23, instead it is correct for the range 
23 < r/kpc < 48. 

Finally, we check how best we can measure j3 in the outer 
most parts, 35 < r/kpc < 84, using the Dll DF and the data in 
hand. This region consists of 762 stars, and we assume a = 4.5 
and 7 = 0.35 and repeat the analysis as in D12. The likelihood 
distributions of model parameters are shown in Figure [9] The 
mass-anisotropy degeneracy is clearly visible here, suggesting 
that it is very difficult to directly measure (3 unless an explicit 
or implicit assumption about the potential is made. 

5. CIRCULAR VELOCITY CURVE OF THE GALAXY 



Here we use the measured values of o>(r) and [3(r) from 
our analysis given in Figure [3] to determine the circular ve- 
locity curve of the Galaxy (v c i rc ) through the Jeans equation 
(Equation [7). Besides anisotropy and radial velocity disper- 
sion information we also need to adopt some density profile 
for the tracer population but not of the spectroscopic sample. 
We adopt a value of a = 2.4 as suggested by re cent works of 
IWafkins et ail (120091) and lDeason et al.l (l2011bl) . for the range 
of distance explored here (r < 25 kpc). The blue dots with 
error bars in Figure \W\ are our estimates of v c h- c using CME 
with ttbin = 1200. The uncertainties on v c j lc were computed 
using a Monte Carlo based scheme from uncertainties in a r 
and (3. The leftmost and rightmost points have comparatively 
larger error bars as compared to intermediate points. For the 
leftmost point the large error bar is mainly due to large error 
in the value of o>. On the other hand, for the right-most point 
the large error bar is mainly due to large error in the values of 
ag and (70. 

In Figure [10] it can be seen that the circular veloc- 
ity profile derived from our analysis (blue points) dis- 
play prominent features. We now check if such fea- 
tures are also observed in other studies using tracers other 

tha n BHB sta rs. For this we over-plot v c ; lc compiled 

by ISofue et all (120091) as black dots, o b tained from sev- 
eral references (iBurton & Gordonl 119781: iBlitz et al.l 11982 ; 
Clemensl fT985t iFich et al.l Il989t iHonma & Sofud Il997bla: 



Honma et al. [ 120071: iDemers & Battinellil 120071) . Further de 
ta ils about the source of each individual point can be found 
in ISofue et ail (l2009|) and references therein. Note, v c ; lc val- 
ues of ISofue et al.l (120091) are computed for (R©, Vlsr) = (8.0 
kpc, 200.0 kms -1 ). Correcting the v irc f° r our adopted values 
of (R Q , Vlsr) = (8.5 kpc, 220.0 kms -1 ) is beyond the scope 
of this work and hence we simply over-plot these published 
values in Figure [10] It can be se en there is a prominent dip 
at 9 kpc in the ISofue et al.l (120091) compiled v c ; lc profile. They 
explain this dip by introducing massless rings on top of a disk 
with exponential surface density. We also find a similar dip 
in our v c i lc profile but at around 10-12 kpc. The slight shift 
in the position of dip could be due t o large width of ou r bins 
and also due to the fact that unlike ISofue et ail (120091) . who 
measure v c i rc in the mid-plane of the Galaxy, we measure v C i rc 
over a spherical shell that intersects with the SDSS footprint. 

We now estimate the mass of the dark matter halo of the 
Milky Way, assuming a three component model of the Galaxy 
consisting of the bulge, the disk, and the halo. The bulge 
is modeled as a Hernquist sp here and the disk is assumed to 
follow an exponential profile (IXue et al.ll2008l). The parame- 
ters for the bulge and disk a re taken from Sofue et all (120091) 
and are kept fixed. Although lSofue et all (120091) use massless 
rings, we here have ignored them since our main aim is to fit 
the dark matter halo. We mo del the dark matter halo using 
the NFW dNavarro et al.ll 19961) density profile. Here we con- 
sider both the halo and the bulge to be spherically symmetric. 
The non-axisymmetric effect due to a bar shaped bulge is ne- 
glected here. Potentials for bulge, halo, and exponential disk 
can be expressed as, 



$diskO) = 
3>bulgeO) = 

$nfw(>") : 



GM disk (l- e -''/ fc ) 



GM bu 



lge 



and 



r+a 

GM vi[ \n(l + rc/R vh ) 
g(c)r 



(10) 
(11) 
(12) 



10 



350 



300 



250 



200 



o 150 



> 



100- 



Exponential disk 
-n Hernquist Bulge 

- NFW Halo 

- Tota 

Observed Average Sofu et al. 2009 




15 20 

r (kpc) 



Figure 10. Circular velocity curve of the Galaxy and their individual components along a galactocentric distance (r). The blue marker represents the value of 
v c ir c obtained in the CME bins in r. Red solid line is our fit of the total potential. Black dotted and dotted-dashed lines are the fixed disk and the bulge circular 
velocity profile for set of ad opted values of masses and scale radii. Dashed line is the fitted NFW profile. Black dots with error bars are the collated Vciic values 
given bv lSofue et al] <2009l) whereas yell ow solid line is the average of the given observed values. 



where M disk = 6.5 x I P 10 M n , b = 3.5 kp c, M bulge = 1.80x 10 10 
M Q , and a = 0.5 kpc (ISofue et al.ll2009h . 

Note, the disk potential as given by Equation (TTOb is spher- 
ically symmetric. It means the disk is considered to be a 
spherical body with exponential surface density fall-off. To 
get an idea on the error that is incurred due to the assump- 
tion of the disk as a spherical body with the mass same 
as the flattened disk we refer reader to iBinney & Tremainel 
(|200l)(Figure 2 .17). Roughly the maximum error in v c j lc is 
13%, which is at a distance about twice of the disk scale 
length. However, at the larger distances along the mid-plane 
the discrepancy is smaller. In the general case of triplanar 
symmetry (elliptic disk), in reality, the disk potential has to 
be the function of both polar and azimuthal coordinates and 
in the special case of axial symmetry (circular disk) it has to 
be the function of sole polar coordinates, in addition to the 
radial coordinates (r). We here use the spherically symmetric 
form for two reasons. Firstly, we make use of the the spher- 
ical form of the Jeans equation given by Equation (|7]i which 
demands a spherical potential. S econdly, it is to ea se the com- 
parison with earlier studies , e.g. JXue et al.l (12008). that adopt 
a similar definition. However, later on we consider a 3D disk 
potential and discuss its consequences on the estimation of 
mass. The function g in NFW potential is given by 



g(c) = ln(l+c)- 



and 



2M vir G 



1/3 



The total potential $(r) of the Galaxy is then simply 



<S>(r) = $bulgeO) + $diskO) + $nfwO) 

the value of Hubble constant, H 



(13) 
= 70.4 

n m = 0.3 dKomatsu et al.llMTTl). and A th = 340 



We adopt 
kms~ 1 Mpc~ 1 . 
dBrvan & Normanlll998l) . 

A NFW halo has two free parameters the mass M v 



and the 



concentration c. Since we do not have enough data points 
spanning a wide range in radius, we avoid fitting both the 
parameters simultaneously. Instead we use the concentration 
mass relation, 



c= 327.3M;;; 1 -, 10 u Mvit/Mq < 10 1 



(14) 



1+c 



as has been reported in N body simulations of dark matter 
halos dMaccio et al.ll2007l) . 

Finally, we can derive the resultant circular velocity 
(vdrc) from the total potential (Equation [T3l by computing 
{rd<&/dr) l l 2 . We fit the obtained theoretical rotation velocity 
curve to our observed values of v c j rc shown by the blue dots 
in Figure [10] and the red line is our best fit circular velocity 
curve. The v crrc profiles for the different components are also 
shown separately. The dashed black line is the corresponding 
best fit NFW halo profile. The best fit value for the fit param- 
eter, M vrr , for our three component baryon and dark matter 
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mass distribution is 0.9^ x 10 12 M Q . 
values of R v ; r and c derived from the best fit value of M vtr are 



The corresponding 



249^3j kpc and 12. (Qj 5 respectively. We estimate the mass of 
the Galaxy within r < 25 kpc to be 2.1 x 10 11 Mq. Assuming 
a functional form for P(vi os /v c i rc ) obtained from simulations 
IXue et all (120081) derive the v c j rc from the v; os of BHB stars. 
The derived v c j rc is then used to estimate the virial mass of the 
dark matter halo. They find M vil - = 0.91^ f & x 10 12 M Q which 
is in good agreement with our result, note uncertainties are 
however, slightly larger in our estimates. Since unlike them 
we do not make any assumption about the functional form of 

P(vios/Vchc)- 

Here we study the effect of chosen disk models for which 
we con sider more realistic three-di mensional potential for the 
disk by iMivamoto & Nagal (1 1975b . which holds for the spe- 
cial case of a circular disk. The expression for this potential 
is given by 



$disk(fl,z) : 



GM d] 



isk 



y / R 2 + (a + Vz 2 + b 2 ) 2 



(15) 



Here again, the disk parameters are obtained from the best 
fit values along the gala ctic mid-plane (z = 0) which repro- 
duces VcircCR ) profile forlSofue et al.l (120091) razor-thin expo- 
nential disk dFreemanlll970l) i.e. b = 0.0 kpc B 

Hence the 

best fit value for the disk parameter a = 2.5 kpc whereas mass 
M disk = 6.5 x 10 10 M Q is taken same as in lSofue etai] ( 120091) . 
Since we assume the disk potential to be three-dimensional 
here, for the purpose of computing the total v c ; lc we con- 
sider the component of force along the radial direction (r) 
only. The bulge and the halo models are kept same as earlier. 
The best fit values for the NFW halo parameters with the up- 
dated disk model is found to be M vil - = 1 .2+^ 4 x 10 12 M Q with 
R v i r = 274^Qkpc. Instead if we consider the total magnitude of 
the force in order to compute the v c ; lc we estimate M v i r (/? v i r = 
269!^ kpc) = 1.1^ x 10 12 M Q . Historically, l arge values of 
a have been assumed wh i le mo deling the disk (IWolfire et al.l 
(1 1 9951) : iBland-Hawthornl (120091) ) . We find that this leads to a 
much more massive dark matter halo- for a = 4.5 kpc we find 
M vir (tf vil - = 299!^ kpc) = 1 .6±^j x 10 12 M Q and for a = 6.5kpc 
we find M vir (7? vir = 321+^kpc) = 1.9^ x 10 12 M . 

There are few things which have insignificant or unexplored 
effects on our mass estimation e.g. we do not take into ac- 
count the mass of the super-massive black hole, which is 
~4x 10 6 M m (ISchodel eTaIll200l lEisenhauer et ai1l200l 
iGillessen et al.l 120091) . and is approximately 1/1000 of the 
mass of the bulge. Moreover, its effect is like that of a point 
mass and can be easily absorbed into the bulge mass. Another 
effect that is not considered is the tidal effect on the primary 
object that has been qualitatively studied to find the mas s ra- 
tio between the Galaxy and M3 1 dBaiesi Piilastrini| | 2009t) and 
depen ds strongly on different impact parameters (lEneev et al.l 
119731) . We find that R vir of the Galaxy is ~ 250 kpc and 
for M 31 it is ~ 260 kpc ( ISeigar et all 120081: iMaiewski etai] 
120071) . Given that the distance between these two galaxies 
is ~ 780 kpc dRibas et al.l l200l iMcConnachie et aD l200l 
iKarachentsev & Kashibadzd 120061) . which is more than the 
double of the sum of their virial radii, we believe the tidal 
effect of M3 1 on overall mass estimate of the Galaxy, if any, 
should be negligible. The tidal effects of LMC and SMC on 
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Figur e 11 . Dashed black line is our best fit model of v'cir e given by red line in 
Figure [Tol Blac k dots are the values from the literatures (Wilkinson &.Evanjj 
1999; Xue et al. 2008; Sofue et al. 2009; Watkins et al. 2010; Gnedin et al. 
l2010; ISamurovic & Lalovic 201 1 ; De ason et alJ2012bl) labeled as W99, X08, 
S09, W10, G10, SI 1, and D12b respectively. To make the plot less obscure 
we do not include similar findings from the literature. For details about the 
similar results refer to the text. 

the Galaxy have not been explored in this paper. 

Note, our mass modeling of the Milky Way, does not make 
any assumption about the value of (3, instead we use the value 
of /3 directly computed from the data. The only assumption 
that we make is that the density of the dark matter halo fol- 
lows an NFW profile. As long as that assumption holds our 
estimates for v c j rc and mass of the Milky Way should also be 
valid in the outer parts r > 25 kpc where we cannot directly 
measure /?. If one wants to directly measure v c j lc in the outer 
parts using only line of sight velocities then one has to make 
an assumption about the underlying /?. Several attempts have 
been made in this regard, with each of them making different 
assumptions about /3 and hence introducing a bias in the esti- 
mated mass. Below we compare these with our prediction for 
the mass of the Galaxy. Dashed line in Figure QT| is our best 
fit model of v c i rc . Since literature sources mostly report mass 
within a certain radius, to facilitate comparison we convert it 
to v C irc using the relation 



M(< r) : 



The disk model with / 



= is also otherwise known as Toomre's model 



■ Kuzmin disk (Toomre 1963; Binnev & Tremaine 2008) 



In Figure QT| the plotted v C irc from different sources span 
a wide range in radii and were computed using different 
types of tracer populations. By fitting a model to the 
kine matics of the satellit e galax ies and the globular clus- 
ters, IWilkinson & Evansl (119991) measure the mass to be 
M(50kpc) ~ 5.4±3g x 10 11 M . This agrees with estimates 
bylKochangj (Qlp) (M (50 kpc) = (4.9 ± 1 . 1 ) x 10 1 1 M Q ) and 
Sakamoto et ail (120031) (M(50 kpc) = 1.8-2.5 x 10 11 M ). 
Watkins et alj (|2010) apply the tracer mass estimator formal- 
ism to 26 satellite galaxies and find that M(300 kpc) = (0.9 ± 
0.3) x 10 12 Mq. Their mass estimate is however prone to the 
systematics introduced by assumed /3, as duly mentioned by 
them. It can be seen in the figure that at r = 100 kpc, depend- 
ing upon the chosen anisotropy, their mass could vary any- 
where between 0.3 x 10 12 M & and 1.4 x 10 12 M Q . Studying 
the hyper-velocity stars within 80 kpc, and assuming (3 = 0.4, 
iGnedin et alj d2010l) estimate M = 6.9^;g x 10 11 M@, which 
is slightly higher th an our estimate. Using BHB stars and 
IWatkins et alj d2010l) tracer formalism ISamurovic & Lalovid 
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Figure 12. Velocity dispersions and anisotropy profiles of BHB stars in sim- 
ulated stellar halo: From top to bottom are q>, ag, aj,, and /3 profiles of the 
1 1 instances of simulated halo taken from Bullock & Johnston 1 2005). The 
thick red lines are the mean profiles of the 1 1 halos. 



( 120111) estimate the mass at r = 85 kpc to be 8.83 ± 0.73 x 10 1 1 
Mq . This is slightly higher than our estimate, probably be- 
cause they assume (3 = 0. With the mixed sample of trac- 
ers (BHB and CN stars) populating t he outer-most halo (r ~ 
50- 150 kpc), iDeason et alj d2012bl) estimated mass of the 
Galaxy to be M(150kpc) = (5-10) x 10 11 M . The varia- 
tion is mainly due to uncertainty on the adopted potential and 
density slopes, and anisotropy. Their range of mass at the 
outer-most halo falls within our estimation. 

6. KINEMATICS OF THE SIMULATED STELLAR HALO 

We now study the kinematic properties of simulated stel- 
lar halos in which the halos are formed purely by accre- 
tion of satellite galaxies. F or this we use the simulations of 
iBullock & Johnston! ( 12005b . In order to construct a synthetic 
sample of B HB stars from these simulations we use the code 
GALAXIA dSharma et al.ll201 lal) . Figure[T2]shows the veloc- 
ity dispersion and anisotropy profiles of 1 1 different ACDM 
halos as a function of galactocentric radius r. The mean of all 
the halos is also shown alongside as thick red line. In gen- 
eral the velocity dispersions fall off with radius. At small r 
the fall is rapid but at large r it is much slower. Asymptoti- 
cally the ratio <7>/V v ir approaches a value of around 0.8. The /3 
rises rapidly from a value of zero in the center to about 0.8 at 
r ~ 10 kpc, and thereafter shows very litt le change. These re - 
sults are in good a greement with results o flAbadi et alj ( 120061) : 
ISales et aT] ( 120071) . who study the stellar halo formed in cos- 
mological hydro-dynamical simulations including star forma- 
tion and feedback. 

Firstly notice that ACDM halos are rarely tangential for any 
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Figure 13. The mean j3 profile of B HB stars in 1 1 simulated ACDM stellar 
halos of Bullock & Johnston 1 20051) . Sh own alongside as d ashed line is the 
best fit analytic function of form given by Cuddeford 1 1991). 
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Figure 14. The fj profile of simulated stellar halos having non-standard ac- 
cretion history. 

given radius. For most of the range of r, j3 i s in general greater 
than 0.5. We fit an analytic function of ICuddefordl (1 1 99 lh 
form given by 



(16) 



to the mean (3 profile of the 1 1 ACDM halos. The best fit 
values of the free parameters were found to be /3q = 0.765 
and ro = 0.00843/? v ir- Figure Q~3] shows that the fit is quite 
good for a wide range of r. The slight mismatch at r < 1 
kpc could be due to issues related to force resolution. In the 
outer parts most of the mass is in bound structures and hence 
is not smoothly distributed. This is probably responsible for 
the non-smooth behavior in (3 in the outer parts. 

Figure [14] presents the beta profiles for halos having non- 
ACDM accretion history, i.e., halos having accretion his- 
tory significantly different from that predicted by the ACDM 
model of galaxy formation. Six halos that we consider are 
with accretion events being dominated by 1) radial orbits 2) 
circular orbits 3) old events 4) recent events 5) higher lumi- 
no sity and 6) low lumino sity and these are from simulations 
of iJohnston et al.l (120081) . Signatures of different accretion 
events can be seen in the (3 profiles. The most significant dif- 
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ference is between the radial and the circular halo, which is 
expected since the orbital properties of the satellites were dif- 
ferent to begin with. It is interesting to note that the circular 
halo is the only one among all the simulated halos that can 
have /? < 0. The old halo almost perfectly follows the mean 
profile that we had derived for the 1 1 ACDM halos and has the 
smoothest profile. This is due to the fact that the stars in this 
halo are completely phase mixed and have no structures of 
any kind. The young halo has very few stars in inner regions 
and shows non smooth behavior due to presence of significant 
amount of structures. The high luminosity halo is also very 
similar to old halo. However, the low-luminosity halo has 
higher /3 for r < 0.01/? v i r . This is most likely due to circular- 
ization of orbits when acted upon by dyna mical friction . Orbi t 
circularization has also been reported by ISales et aT] ( 12007b 
in their simulations. Note, the effect of dynamical friction 
is strongest for high luminosity events and weakest for low 
luminosity events. Moreover, satellites when acted upon by 
dynamical friction loose energy and move towards the inner 
regions of the halo. This partly explains as to why the high 
luminosity halo has low /3 in the inner regions as compared to 
the lo w lu minosity halo. 

In 33.2l we measured /3 until r = 23 kpc. Beyond this we can 
measure a r out to r = 56 kpc, but not a$ and ag. By using the 
circular velocity curve that we derived in we can predict j3 
beyond r > 23 kpc making use of the Jeans equation (Equa- 
tion|7]). To proceed we need to make an assumption about the 
density slope (a); beyond r > 27 kpc it has been shown that 
a is around 4.5. Assuming this, the predicted j3 is plotted in 
Figure[l5] It can be seen that there is a slight jump in the value 
of f3 passing from r = 23 kpc to r = 27 kpc and beyond this the 
value of j3 is around 0.4. The sudden jump in j3 occurs via 
the Jeans equation (Equation[7]l, due to the discontinuity in a 
(=-dlnp/dlnr). Note, an assumption of steeper density slope 
would increase the (3 and vice versa. The red line in figure 
is the anisotropy profile fitted to the simulated ACDM halos 
given by Equation ( fTol l from $6] It can be seen that accretion 
based models cannot explain the dip that is in observations, 
specially the profile in the region 12 < r/kpc < 23. However, 
outside this region the simulations are roughly in agreement 
with observations. For r < 12 kpc the observations match the 
value of around j3 = 0.5 seen in simulations. For r > 23 kpc 
the overall value of (3 in observations is slightly low but the 
profile is flat as in simulations. The outer halo at r = 56 kpc is 
radial with (3 = 0.55. 

7. CONCLUSION AND DISCUSSION 

We study the kinematics of ~ 4500 BHB stars to obtain 
the velocity dispersion profiles along three orthogonal axes in 
spherical polar coordinates using the gaussian velocity ellip- 
soidal (GVE) DF. GVE as an estimator of the velocity disper- 
sion has the advantage that no assumptions about potential or 
density are needed a priori. From the estimated velocity dis- 
persion profiles using maximum likelihood analysis, we also 
derive the anisotropy profile of the Galactic halo and compare 
it to the simulated ACDM halos. Finally, using radial velocity 
dispersion profile, anisotropy profile and density power law 
we constrain the mass of the Milky Way Galaxy using the 
Jeans equation. 

We measure the a r profile of the halo out to r ~ 60 kpc. 
At large distance (d 3> Rq), o> can be approximated by oi os . 
Thus in outskirt, eri os profile given bv lXueetalJlHOOl con- 
verges to our oy profile. At r ~ 60 kpc a r attains 100 kms -1 . 
However, in the inner halo the approximation (o> m a\ os ) is 



invalid and we find that the deviation of eri os from a r is as 
high as ~ 40 kms" 1 . We obtain a o>(r) profile with plateau in 
the inner halo, a sudden fall at r ~ 15 kpc and a gradual de- 
cline outwards. Quali t atively , similar profile is also found by 
ISommer-Larsen et al.l (119971) . However, our a r profile sharply 
falls at r ~ 15 kpc whereas they find a gentle transition. The 
is probably due to the fact that they make an assumption that 
VdrcM is constant which we find is not completely true. 

Next we estimate the tangential velocity dispersions, ag and 
<70. Using these estimates we are able to measure the j3(r) till 
r = 25 kpc. Astonishingly, we discover a dip in the j3 profile at 
r = 17 kpc, where j3 m — 1.2. We find that the inner halo (r < 12 
kpc) is radial with f3 ±2 0.5. This result o f radially biased in- 
ner halo concurs with the recent results bv lSmith et al.ll2009bfc 
iBond et alJ20Tol using the proper motions. Beyond the switch 
over point in the range 18 < r/kpc < 25 the anisotropy rises 
slightly and becomes isotropic to mild radial. We also verify 
the result using an alternative DF, namely the Dl 1 DF. A small 
systematic in the f3(r) profiles is seen from these two models 
which is mainly due to the assumption about density slope 
(a) that needs to be made a priori for Dl 1 DF. We check for 
the contribution of the halo substructures, namely Virgo Over- 
density and Sagittarius stellar stream, and find that they have 
little effect in the anisotropy profile. The effects of Vlsr and 
Rq upon our velocity dispersions and anisotropy estimates are 
also found to be negligible. 

Dl 1 study the BHB stars in the radial bin (10 < r/kpc < 25) 
and find the halo to be tangential. After re-analyzing the Dl 1 
sample in this bin, we find that this is mainly because of the 
choice of their bin size that encompasses the transition re- 
gion (13 < r/kpc < 17) where we detected a dip in the (3(r). 
Possibly it could be also because of the degeneracy between 
potential and anisotropy in their model. However, in their 
recent work D 12 break the degeneracy among their model pa- 
rameters and measure (3 = 0.4±Q2 5 with a = 4.6 in the region 
16 < r/kpc < 48. Within the range of uncertainty our value 
for f3 (— 0.14^Qg?) using GVE model agrees with them. We 
find that D12 value of j3 « 0.4 in this bin, although derive 
from a sample which is dominated by stars within r < 27 kpc, 
is not appropriate for the range 18 < r/kpc < 23, instead it is 
appropriate for the range 23 < r/kpc < 48. Finally, we check 
how well we can estimate the mass and anisotropy together 
using Dl 1 DF in the outer- most region (35 < r/kpc < 84) We 
find that due to the lack of tangential information, a degener- 
acy between mass-anisotropy cannot be broken. 

Substituting the estimates of o>(r) and j3(r) in the Jeans 
equation, we then calculate the circular velocity profile of 
the Galaxy (v c i rc (r)). We det ect the dip in the v c ; rc pro- 
file at 10-12 kpc, also seen bv lSofue etaTl (120091) at 9 kpc, 
which is attributed to the massless ring as a perturbation 
to the disk. Finally, we fit the three component (exponen- 
tial disk, Hernquist bulge, and NFW halo) galaxy model to 
the observed v c jr C profile in order to obtain the mass dis- 
tribution of the Galaxy. From our best fit model, we cal- 
culate Mvii of the halo to be 0.9^ x 10 12 M with R vir = 
249^31 kpc and concentration parameter, c = 12. O^. The 
mass of the Galaxy, within the extent we are able to con- 
strain all the three components of velocity dispersions, is 
estimated to be M(r < 25 kpc) = 2.1 x 10 11 M Q . Our es- 
timate for M v ; r is in good agreement with the most of the 
recent e stimates namely bv lWatkins et alj d2010l): iKochanekl 
(119961): IWilkinson & Evans! (119991): ISakamoto et al.l (120031): 



ISofue etalJ (12009b : IWatkins et alj (120101) : iDeason et al] 
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Figure 15. Black solid line is the observed profile in the outer halo (assumed a = 4.5) whereas the red solid line is the ft profile of the simulated halo. The 
blue dots with error bars are the observed values of f3 from section [3] The sudden jump in profile passing from r = 23 kpc to r = 27 kpc could be a spurious 
effect due to our assumption of the broken power law with break at r = 27 kpc. 

Alternatively, the shift in the properties seen in the ob- 
served profiles could possibly be an indication of a com- 
plex multi-component halo. Recently there have been series 
of works advocating multi-component hal o. The studies of 
the cal ibration stars by lCarollo et al.l (120071) and lCarollo et al.l 
(l2010l) from the S PSS survey and the follow-up studies by 
iBeers et all (120121) have shown that the halo has at-least two 
distinct components. They associated the inner-halo to be 
formed in-situ where as the outer halo are con sidered to be 
forme d by accretion. ICarollo et al.l (120101) and Ide Jong et al.l 
(120101) have found that the population fraction inversion point 
between the inner an d outer halo lies between ~ 15-20 kpc. 
iKinman etafl ( 120121) studied the BHB and RR Lyrae stars 
towards the galactic-anticenter and North Galactic Pole and 
found that the retrograde component of the halo dominates for 
r > 12.5 kpc. It seems that this transition between the inner 
to the outer halo is recorded in the f3 of the halo as well. Ad- 
ditionally, duality in the formation history of the halo has also 
been seen in the recen t smooth particle hy d rodynamics and 
N-bod y simulation s by IZolotov et alJ (120091) ; iMcCarthv et al.l 
(120121) ; iFont et al.l (12011b . iKalirail (120 1 2b recently attribute 
the age difference of 2 billion years in the halo components to 
in-situ and accretion. On the contrar v.|Schonrich et al.l (1201 lb 
reanalyzed the calibration stars from lCarollo et al.l (12010b and 
find no reliable evidence of the existence of outer retrograde 
halo. 

In a nutshell, the stellar halo is a test bed to understand the 
formation history of the galaxy. Even in this era where we 
have access to huge volume of spectroscopic and photomet- 
ric data, the crucial physical quantities like velocity disper- 
sions and anisotropy are not completely understood due to the 
lack of proper motions. With the advent of data inflowing in 
the coming decades through the magn ificent next gener ation 
of spectroscopic survey like LEGUE dDeng et alj|2012b and 
specially, unprecede nted proper moti ons from an astrometric 
mission like GAIA (lPerrymanll2002b will help to put strong 



d2012bb as demonstrated in Figure ITT] In t heir st udies of 
same population of stars (BHB), IXue et all (12008b also fit 
a three component galaxy model and calculates M v j r to be 
0.91^0 ig x 10 12 M Q . Our result forM v i r is in very good agree- 
ment with their estimate but our uncertainty is slightly larger. 
Note that they make an assumption about the (v irc/vi O s) with 
radius from the simulations whereas we do not make any such 
assumption. Additionally, we also consider a more realis- 
tic three-dimensional disk model which is found to predict 
slightly higher Galactic mass M vir = 1.2!°;|x 10 12 M Q with 
i? v ir = 274+3Q kpc for flattening constant a = 2.5. 

In the end, we used the measured quantities o> and v C ir C to 
extend the /3 profile beyond r ~ 25 kpc up to the distance 
where cr,-(r) can be confidently measured (r ~ 60 kpc). The 
only assumption that we make here is about the density profile 
w hich we choose to be a = 4.5 in consent to the recent results 
by IWatkins et all (12009b : iDeason et al.1 (1201 lbb . We find that 
the outer halo is radial and attains j3 = 0.55 at r ~ 60 kpc. 

We also compare our result with si mulated stellar halo 
which are formed purely by accretion dBullock & Johnston! 
12005b . This simulated halo is found to be in rough agree- 
ment to the observed halo in the inner region r < 12 kpc. It is 
seen that in none of the instances of simulations the j3 profiles 
obtained could predict tangential halo at any distance and thus 
fails to explain a dip seen at r = 17 kpc in observed /3 profile. 
In contrast, in the outer region (r > 25 kpc) simulations and 
observations both agree in overall sense of the anisotropy and 
predict a flat anisotropy profile. 

In all the observed quantities <r r ,<rg,a^ and /3 we see a dra- 
matic shift in properties at r ~ 17 kpc. We noticed that these 
undulations in the profiles are translated into our v c i, c estima- 
tion resulting a varying v c i lc profile. It could be true other way 
around, in a sense that the non-monotonic trends seen in all of 
our kinematic profiles could be due to the presence of so far 
unaccounted features in the Milky Way potential. 
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constrains on these fundamental quantities. Additionally, to 
see a bigger picture, confirming the results with the different 
stellar types or an alternative tracers is also crucial. More- 
over, exploring the southern sky is equally important to com- 
plete the picture, for which up-coming spectroscopic survey 
like GALAH0 will also play an important role. 
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APPENDIX 

A. BINNING AND EFFECT OF THE BIN SIZE 

Here we investigate an effect of «bin to our analysis. Figure[l6]shows the velocity dispersion profiles and the anisotropy profile 
for the same sample of stars but with different particles in each bin. We see that with the decrease in the number of stars in each 
bin the uncertainties in the result increases. But the overall trend remains unaltered. 

B . EFFECT OF V LSR AND R 

In literatures there are varied claims about the value of Vlsr ranging from 1 84-270 kms -1 (lOlling & Merrifieldl 119981 : 
iMendez et all 1 9991: iReid et al.ll2009UBoyv et al.ll2009l:lKoposov et al.ll2010l). Similar l y, the value of R is also disputable within 
8-8.5 kpc (lReiJl993HGhez et alJl2008t iGiUessen et al.ll2009l) . iMcMiUan & Binneyl dSoToh found that the ratio v LS r/R© can be 
better constrained than each of them alone and should range between 29.9-31.6 kms _1 kpc _1 . Distressingly, there is still no 
consensus upon the values of (vlsr, Rq)- To study the effect of chosen values of (vlsr, Rq) upon our estimates of dispersion 
profiles we repeat the same analysis done to obtain the black diamond points in §|3] (Figure |3j. In Figure[l7]we show the results 
for different values of (vlsr, Rq)- Here again the black diamond markers are obtained for a case (vlsr, Rq) = (220.0 kms -1 , 8.5 
kpc) and is thus a replica of diamond points from Figure |3]put again for the ease of comparison. The red and black markers in 
the figure show the effect of chosen R Q values upon our estimates whereas the red, blue and cyan markers demonstrate the effect 
of chosen Vlsr- Within the range of (vlsr, Rq) investigated, the figure clearly demonstrates a negligible effect of them upon our 
estimates, given the range of uncertainties. 



3 |http : / /www . aao ■ gov ■ au / HERMES /GALAH/ Home ■ html| 
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Figure 16. Effect of the number of particles in each CME bin. Red, green and blue points are the estimates with 750, 1000 and 1200 number of stars in each bin 
respectively. From top to bottom are the ay, erg, cr^ and anisotropy profile. 
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Figure 17. Velocity dispersion and anisotropy profiles for different combinations of vlsr and Rq . 
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